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Abstract 

We study a phase transition in a non-equilibrium model first in- 
troduced in 0, using the Yang-Lee description of equilibrium phase 
transitions in terms of both canonical and grand canonical parti- 
tion function zeros. The model consists of two different classes of 
particles hopping in opposite directions on a ring. On the complex 
plane of the diffusion rate we find two regions of analyticity for the 
canonical partition function of this model which can be identified by 
two different phases. The exact expressions for both distribution of 
the canonical partition function zeros and their density are obtained 
in the thermodynamic limit. The fact that the model undergoes a 
second-order phase transition at the critical point is confirmed. We 
have also obtained the grand canonical partition function zeros of 
our model numerically. The similarities between the phase transi- 
tion in this model and the Bose-Einstein condensation has also been 
studied. 

Key words: Matrix Product Ansatz, Asymmetric Exclusion Pro- 
cess, Yang-Lee theory 
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1 Introduction 

One of the most important activities in the field of equilibrium statistical 
physics is the study of the phase transitions; nevertheless, there is a gen- 
eral framework for the statistical description of the equilibrium systems 
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and also different approaches for studying their equihbrium phase transi- 
tions. One of these theories was proposed by Yang and Lee in 1952 
The Yang-Lee theory of equilibrium phase transitions is based on the zeros 
of the partition function. It is known that the zeros of the grand canonical 
partition function of finite systems Z{z) (in which z is fugacity) are gener- 
ally complex or negative if real. In the thermodynamic limit roots might 
move down (or up) and touch the positive real axis. When this happens, a 
phase transition occurs because the system can have different behaviors for 
z < Zq and z > zq, where zq is the value of the root on the real axis. For 
example the pressure P = kBT\miv~>ooii^/V) logZ) will be non- analytic 
and the density p = {d /d log z){P/kBT) will be discontinuous at the tran- 
sition point which predict a first-order phase transition. Similarly, one can 
investigate the zeros of the canonical partition function as a function of 
complex-temperature and find the same transition points 0]. By calculat- 
ing the line of zeros and also their density in the thermodynamic limit one 
can find the transition point and also the order of transition exactly. 
Recently much attention has been focused on one-dimensional out of equi- 
librium systems because of their interesting properties such as first-order 
phase transitions and spontaneous symmetry breaking jDIS- However, in 
contrast to the equilibrium systems many powerful concepts are missing in 
this context. For example, the applicability of the Yang-Lee theory to the 
non-equilibrium systems such as one-dimensional driven diffusive models 
is a quite non-trivial question and yet without answer. People have tried 
to apply the Yang-Lee theory to describe phase transitions in these models 
IHIIZI- It seems that one can define similar quantities such as a grand canon- 
ical partition function and then apply this theory to the out of equilibrium 
systems without any problem. In this paper we will apply the Yang-Lee 
theory to an exactly solvable one-dimensional non-equilibrium model and 
investigate its phase transitions. This model has already been solved in 
our previous work and exact results are available Here we are going 
to compare our previous results with those obtained from application of 
the Yang-Lee theory; however, we should mention some of the differences 
between our work in this paper and what other people have done so far. In 
|5] the author studies a particle-conserving driven diffusive model consists 
of two different classes of particles with finite densities ^ ^ and applies 
the Yang-Lee theory to it. He aimed to locate a phase transition induced 
by varying the density of particles. By introducing a grand canonical par- 
tition function as a function of a fugacity-like quantity a first-order phase 
transition is located by studying the numerically obtained zeros of this 
function. In another paper (21 the authors investigate the phase transitions 
in the asymmetric simple exclusion process (ASEP) with open boundaries 

^This model is known as AHR model in related literature. 
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[TUl [TT] by studying the zeros of the partition function of the model. In 
their approach the boundary rates were generahzed to the complex plane. 
They obtained the distribution of zeros near the transition point and also 
their density near the real axis using the similarities with electrostatic the- 
ory. They could also calculate the line of zeros analytically. 
The solvability of our model allows us calculate its canonical partition func- 
tion exactly. This was done in our previous paper |S] and will be reviewed 
in the second section. Apart from numerical estimates for the zeros of the 
canonical partition function as a function of reaction rates, we have calcu- 
lated the line of zeros and also their density (which determines the order of 
the phase transition) analytically using the equilibrium statistical physics 
toolbox. Having the exact analytical results we can discuss the possibility 
of phase transition and also obtain its order. Next we will define a grand 
canonical partition function and study its zeros in the complex-fugacity 
plane. The properties of this function reveal the similarities between the 
phase transition in our model and that of a Bose gas. In the last sec- 
tion we will summarize our results and generalize our approach to other 
non-equilibrium models. 



2 The Model 

In [S] we introduced a one-dimensional exclusion model consists of two dif- 
ferent classes of particles (we call them positive and negative particles here- 
after) which occupy the sites of a chain of length L with periodic boundary 
condition. Each site of the chain is cither empty or occupied by a negative 
or by a positive particle. The positive (negative) particles hop to their 
immediate right (left) sites with unit rate provided that the target sites are 
empty. Adjacent particles with different charge type might exchange their 
positions with asymmetric rates 1 and q. Specifically, the interaction rules 
are 

+ — > + with the rate 1 
— — > — with the rate 1 , , 

+ — — > — + with the rate q ^ ' 

— + — > + — with the rate 1. 

Assuming that there are M positive particles and only one negative particle 
on the chain ^ we showed that the steady state weights P{C) of the model 
can be obtained exactly using a so-called Matrix Product formalism 
One can define the sum of these steady state weights as a quantity which 
plays the role analogous to the canonical partition function in equilibrium 

^In this case our model is a special case of AHR model |H| in which the number of 
the positive and negative particles on the ring are equal. 
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statistical physics 

c 

It turns out that Z has a closed form in terms of the transition rate q, the 
number of the positive particles M and the length of the chain L 

i=0 ^ 

in which Cf — ■i^.'^ ...i is the binomial coefficient. In the thermodynamic 



limit 



M 

L , M — > oo with p = — being fixed (3) 



using the steepest decent method it can be shown that 



for q <2p ZL,M{q) ^ (^_2 



/ 2\L-1 

(4) 



for q>2p Zl.m{<D ^ (1 - + )gf ■ 

The existence of two different phases is apparent. In the same reference 
we have shown that the density profile of the positive particles has an 
exponential behavior for q < 2p with a correlation length ^ = |ln^|^^ 
which diverges as q approaches its critical value qc — 2p, while it is an 
error function for q > 2p. This proves the existence of a second-order 
phase transition from a power-law phase to a jammed phase. Using the 
Matrix Product formalism one can also calculate the speed of the different 
species of particles on the ring in the thermodynamic limit. Both speeds 
are linearly increasing functions of q ior q < 2p. However, for q > 2p the 
speed of the positive particles is a constant equal to 1 — p while the speed 
of the negative particle is a complicated increasing function of q. In the 
following section we will apply the Yang-Lee theory of equilibrium phase 
transitions to our model. As we will see, it will not only recover all the 
above mentioned results but also shed more light on the unknown aspects 
of our problem. 



3 The Partition Function Zeros 

Let us study the phase transition of the our model using the Yang-Lee 
theory. We consider the zeros of the canonical partition function of the 
model Z^^Miq) given by ||2J) in the complex-g plane at fixed L and M . In 
Figure (1) we have plotted the numerical estimates of these zeros for a 
chain of length L = 150 and three different values of p. As can be seen 
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Figure 1: The numerical estimates for the roots of ZL.Atiq) in the complex- 
q plane for different values of density p = 0.6, 0.8, 0.9 and L = 150. It is 
seen that the roots accumulate to the positive real q axis at 1.2, 1.6 and 
1.8 respectively. 



the zeros accumulate slowly to the real q axis at a critical value qc — 2p. 
As we will see later this accumulation takes place at an angle ^ which is 
the reminiscent of a second-order phase transition jj. In what follows we 
will try to find the line of the canonical partition function zeros and also 
their density near the positive real-g axis using the equilibrium statistical 
physics tools. We will see that the equilibrium-type calculations give the 
same results obtained from the numerical estimates. 
It has been shown that the line of zeros can be obtained from ^ 

Re gi = Re g2- (5) 

In equilibrium statistical physics g is the extensive part of free energy and is 
generally a function of the temperature (here q) and the density of particles 
p. The indexes 1 and 2 show the values of the function g in the right and 
the left hand side of the critical point. Here we define this function as 

.9(9,P)= lim ]-\nZL,M{q) (6) 

where lim ■ • ■ is in fact the thermodynamic limit given by Q . Using the 
asymptotic behavior of Z^^Miq) given by (0J we can calculate g{q,p) in 
each phase. After further calculations we obtain 

for q<2p g{q,p) ^ In ^^^^^^ 
for q>2p g{q, p) = In pp(i„p)i-p . 
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By substituting Q in Q we find the following equation for the line of 
zeros in the complex-^ plane 

in which x = Re{q) and y = Im{q). In Figure (2) we have plotted both 
((SJ and the numerical estimates of zeros given in the Figure (1) for three 
values of p. As can be seen the curves lie on the numerical estimates 





Figure 2: Plot of the equation © (solid lines) for three values of density of 
the positive particles p = 0.6 (inner), p = 0.8 (center) and p = 0.9 (outer). 
They cross the positive x axis at a; = 2p. The dotted lines belong to the 
numerical estimates of the canonical partition function zeros for the same 
values of densities and are taken from the Figure (1). 

of the canonical partition function zeros. The small difference belongs 
to the fact that the numerical estimates have not been calculated in real 
thermodynamic limit. The curves also cross the positive x axis at Xc = 2/5 
which is the transition point as we had mentioned above. In the equilibrium 
Yang-Lee theory it is well known that the density of zeros on the real 
positive axis is zero at a second-order phase transition. In order to obtain 
the density of zeros p in this region (on the positive real x axis and near 
the critical point) we use the following equation first introduced in ^ 

d 

2ttp{s, p) ^ —Im{gi - gs) (9) 

in which s is the arc length of the line of the zeros which is zero at the 
critical point and increases along with the line of zeros in the positive y 
direction. The values of gi and 172 are given by Q. In order to calculate 
p{s, p) first we find an equation for the line of zeros which is valid for small 
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Figure 3: Plot of the numerical estimates of the zeros of the grand canonical 
partition function H12|l for q = 0.5, 1.0 and 4.0. The length of the chain is 
L = 50. 



y's in the vicinity of the critical point Xc ~ 2/3. Using ((SJ it can be shown 
that this is actually a line 

y = 2p - a; for |a; - 2p| -C 1 , < y < 1. (10) 

The equation H10|l confirms that the accumulation of zeros in the vicinity 
of the real-g axis takes place at angle j. Using (|5J and H10() we find 
for p 7^ 1 

dy d 2 

2My,P) = - a^^-an ^ ^yy^^ „ ^ _ oc (11) 

Therefore, as y approaches to zero the density of zeros becomes zero as 
we expect for a second-order phase transition. Comparing the numerical 
data given in Figure (1) with the results obtained from the application of 
the equilibrium statistical physics tools shows good agreement between the 
two approaches; therefore, we it is reasonable to believe that the analytical 
approach presented here gives the exact result. 

It is also interesting to look for the zeros of the grand canonical partition 
function of our model. In this case we can investigate the similarities 
between the phase transition in our model with that of a Bose gas. As 
opposed to the equilibrium statistical physics, the definition of a grand 
canonical ensemble for the steady state of a non-equilibrium system is not 
unique. We adopt the following definition first used in [T^ 

Zi(g,z)= ^ z*^Zi,M(g) (12) 

A/=0 
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in which z is the fugacity of the positive particles and Zl^m{<i) is given 
by 101 ■ It is seen that the grand canonical partition function H12I) is a 
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Figure 4: Plot of the density of the positive particles H14() as a function of 
fugacity z for L — 50, 100 and 200 at g = 0.3. In the thermodynamic limit 
there is a finite discontinuity in the density of particles. 

polynomial of degree i — 1 in z so, in the complex-z plane it has L — 1 
zeros Zi and can be written as 

L-l 



ZL{q,z)^\[{z 



i=l 



Using (PI one can easily calculate the grand canonical partition function 
(O explicitly 



Z,iq,z) = i^—^)A^-^ 
q A q 



(1 + ^) 



L-l 



(9-3) 



+ 



L-l 



(13) 



The fugacity of the positive particles in H12|l has to be fixed by density of 
them 

p{z)^^^^\nZL{q,z). (14) 

Let us examine the zeros of the grand canonical partition function given by 
(|13ll in the complex-z plane. Since the grand canonical partition function is 
a real polynomial with positive coefficients (they are sum over probabilities) 
the zeros come in complex conjugate pairs and the real roots are negative. 
In Figure (3) we have plotted the numerical estimates of the zeros of p3fl 
in complex-z plane for three values of q. As long as q < 2 the roots lie 
on a vertical elliptic. It also appears that the zeros approach the real z 
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axis at an angle ^ . It is a sign of a first-order phase transition. For q> 2 
all of the roots have negative real parts; therefore, we do not expect any 
phase transition to take place. For q — 2, we should take the limit of (O 
since it is undefined at this point. In the following we will investigate our 
model in order to see its similarities with the Bose-Einstein condensation. 
Using H13|) one can calculate H14|) in each phase in the thermodynamic limit 
explicitly. It turns out that 



^ for z < -2- 
'^^ ' I 1 for z > 



(15) 



In Figure (4) we have plotted H14|l for different values oi L. As can be 
seen the L dependence of the density suggests that in the thermodynamic 
limit, p increases with z and then at a specific point zq = 5^ ^ 
finite discontinuity. At this point, z does not fix the density anymore and 
the system undergoes a first-order phase transition. The finite jump in the 
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Figure 5: The pressure P as a function of p^^ for q = 0.2 and L = 15000. 

density is related to the finite density of roots of (|13|l . /i(z), at the real z 
axis 

1 - p(^o) 



2hzq 



(16) 



The critical fugacity zq can also be obtained by extrapolating the real 
part of the nearest root"^ to the real positive z axis for L — > cxa. In the 
Bose-Einstein condensation the density of the particles has such a behavior 
where the conservation of the number of particles is broken ^3] . Another 
interesting quantity is the pressure which can be defined analogously to 

■^In fact there are two of them since the roots appear in complex conjugate pairs. 
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equilibrium physics 

P{z) = jln ZL{q,z). (17) 

However, P is not the physical pressure of the particles in the current 
context. The particles pressure in each phase in the thermodynamic limit 
can be calculated using l(T^ 

r ln(l + z) for z < 
V ' I in(^) for z > ^ ' 

Now by using |[T5|l and ((T5|) it can easily be verified that in the power- 
law phase (p < |) the particles pressure P as a function of density has 
the form P{p) = ln{Y^) while in the jammed phase {p > |), as can be 
seen from (11511 . the density- fugacity relation 1)14(1 breaks down and results 
in /5 = 1^; however, the particles pressure remains constant in this phase 
P{p) — ln(2Z^)- In Figure (5) we have plotted ((17|l as a function of the 
inverse of density of particles p~^ for L — 15000 and q ~ 0.2. The critical 
density in this case is p = 0.1. As can be seen for ^ > 10 the pressure 
decreases as the density gets smaller; however, for i < 10 it is nearly 
constant. This isotherm (here q instead of temperature T) is similar to the 
isotherm of the free Bose gas when the Bose-Einstein condensation takes 
place. One can also look at the compressibility k which is defined as 

K = LC (19) 

in which 

J{p^)-{pf 

C = — (20) 

P 

and the fluctuation of the density of the positive particles can be obtained 
using 

(p')-(p)' = f|-(f|-ln^L(g,^)). (21) 
L oz L oz 

Figure (6) shows C as a function of the chain length L ior p — 0.6 and 
two values of q above and bellow of the transition point. The transition 
point in this case occurs at qc = 1.2. It can be seen that C is a decreasing 
function of L in the power-law phase {q > qc) while it is nearly constant 
in the jammed phase [q < qc). Using (|13|l and l|21|) it can be verified that 
in the thermodynamic limit ^ drops down as in the power-law phase 
while in the jammed phase it remains constant which gives a divergent 
compressibility. 

^Physically, this is related to the existence of a shock in this phase. 
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Figure 6: Plot of 12()|l as a function of L for p — 0.6 and two values of q. 

4 Concluding Remarks 

In this paper we have used the Yang-Lee description of equilibrium phase 
transitions in terms of the zeros of both canonical and grand canonical par- 
tition function to study a non-equilibrium phase transition in a driven dif- 
fusive system. By studying the canonical partition function zeros a second- 
order phase transition was predicted. This is in quite close agreement with 
our previous results in 5 . The line of zeros and also their density on the 
real axis near the transition point were obtained exactly. By introducing 
the grand canonical partition function of our model (|12() the similarities 
between the phase transition in our model and the one seen in the equilib- 
rium Bose gas were elegantly observed. 

The approach that we used in this paper can also be applied to other mod- 
els. In a similar model consists of a group of first class particles hopping 
behind a slow (or a second class) particle on a closed chain of length L, 
with the following interaction rules 



10 — > 1 with the rate 1 
2 — > 2 with the rate a. 



(22) 



it is shown that by choosing the right reaction rates the probability dis- 
tribution for the stationary state can be mapped to the one obtained for 
an ideal Bose gas 9^. The thermodynamic limit of the canonical partition 
function of this model is obtained in 9 



for a < 1 — p Zl^m{ol) 
for a> I - p ZL,M{a) 



ap' 



OL — 1 L, 



(23) 
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in which p = is the density of the first class particles. Now using © 
and ^ wc find the line of zeros in complex-a plane 

iC-^r + i^-miT^)' + {j^rr-" = i (24) 
p p 1 - p 1 - p 

where x = Re(a) and y = Im(a). This function crosses the real positive a 
axis at Qfc = 1 — jO at an angle ^. It can be seen that the thermodynamic 
limit of the canonical partition function of our model is quite similar 
to (|23|) . Using (|3J) one can also obtain the line of zeros for the ASEP with 
open boundaries exactly. For this model with two parameters a and [3 as 
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Figure 7: Plot of the equation (^5)1 for /5 = -i and /3 = 1. 



injection and extraction rates of particles, the line of zeros in the complex-a 
plane (and fi being fixed) is 

{-x^ +x + y^f + {y-2xyf = .f (25) 

in which we have defined x = Re{a), y = Im{a) and the current of particles 




3 if/3>^ 
/3(l-/3) if/3<i. 



The investigation of H25(l confirms the existence of two phase transitions 
in the system. In Figure (7) we have plotted the line of zeros (|25ll for two 
values of (3. It can be seen that (|25|l gives exactly the same results obtained 
in [7] from the numerical estimates of the canonical partition function zeros. 
In |7] the line of zeros is stated to be |a(l — a)| = J. By putting a = x + iy 
one finds H25|l for the line of zeros. 

Apart from the above mentioned models one can show that our approach 
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can be applied to many other non-equilibrium systems like those introduced 
in Work in this direction is in progress |15|. 
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